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Abstract: Borigi-Creutz (BC) model describing the dynamics of light quarks in lattice 
QCD has been shown to be intimately linked to the four dimensional extension of 2D 
graphene refereed below to as four dimensional graphene (4D- graphene). Borrowing ideas 
from the field theory description of the usual 2D graphene, we study in this paper the 
anomalous quantum Hall effect (AQHE) of the BC fermions in presence of a constant 
background field strength with a special focuss on the case = Be^i + Eexi^v 
with B and £ two real moduli and det T^v = B 2 x£ 2 . First, we revisit the anomalous 2D 
graphene by using QFT method. Then, we consider the AQHE of BC fermions for both 
regular det ^ and singular detF^ u = cases. We show, amongst others, that the 
exact solutions of the BC fermions coupled to constant J-^ v have a 5D interpretation; and 
the filling factor ubc of the BC model coupled to constant is given by 24 ^ 2JV+1 ^ 2M+1 ^ 
with N, M positive integers. Others features, such as T^^ D / and the extension of 
the obtained results to the lattice fermions like Karsten-Wilzeck (KW) fermions and naive 
ones, are also discussed. 
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1. Introduction 



Few years ago M. Creutz made a remarkable observation about links between the physics 
of 2D graphene and four dimensional lattice chromodynamics (QCD) [l], [|]. The two Dirac 
valleys of 2D graphene U-f| have an analogue in the J^D extension of Creutz; and are 
interpreted as the up and down quarks of the physics of QCD |], [?], || ^, 10]. This corre- 



spondence opened a window on applying constructions of 2D graphene modelings to J^D 
lattice QCD formulated on the hyperdiamond lattice Q-|1S]. 



In this paper, we contribute to this matter by looking for the generalization of some specific 
properties of graphene to the case of ^D lattice QCD with a special focus on Borigi-Creutz 



(BC) fermions ]|16|-|21]; in particular the aspect regarding the anomalous quantum Hall 



effect (AQHE) ||22(| -[31|. Recall that QHE in higher dimensions has been first considered 



in |27|] for the case of the 4-sphere and has been further developed in [28, 2S, 30] and refs 



therein for other higher D- manifolds. These models are non relativistic systems general- 



izing the well known 2D Hall fluid [32, ^] where, due to disorder, the conductivity has 
plateaus resulting from the non uniformity of the spatial potential. Recall also that the 
quarks u a (x) and d a (x) are light particles that play a central role in lattice QCD simula- 
tions; they are described by 4D Dirac spinors ^> a (x) with integral color quantum numbers 
as required by the SU C (3) non abelian gauge symmetry; but fractional electromagnetic 
charges respectively equal to Q u = ^ an d Qd = — §• Interactions between quarks are 
mediated by gauge fields; in particular by = A a ™Q em + Aj/Ti valued in the adjoint 
representation of the gauge symmetry G = U em (1) x SU C (3) with the 8 hermitian 3x3 



matrices 7} standing for the basis generators of SU C (3). In the QCD regime p7\ , p8| , 39] 



where the SU C (3) x U em (1) gauge interactions can be approximated by constant back- 
ground fields F^ v = (F^u , ) , one is left with a physics quite similar to the one 
of the AQHE of the delocalized electrons of 2D graphene. Therefore, one expects the light 
quark's dynamics to show as well an anomalous quantum Hall type phenomenon 1 in the 
presence of a constant T^ v . This tensor reads generally in terms of the gauge potential 
as d^A^ — dyA^ + but for explicit computations, we will mainly focus on the 
abelian part of the gauge symmetry G. 

Because of the uniformity condition -§^^ u = 0, and also by restricting J^ v to take values 
in the Cartan subalgebra U% (1) x U em (1) C G, the abelian part of the gauge potential 
A 1 ^ (x) that obeys the Lorentz condition d^A^ 1 = reads, up to irrelevant numbers, as 
follows \T 1 jliV x v . Notice by the way that unlike 2D graphene, the underlying space time of 
4D lattice QCD has four real euclidian dimensions leading to several possibilities for the 
allowed directions of J-^. This richness, which has also a physical interpretation in terms 
of interactions (see section 4), may be fixed by looking for configurations that permit the 
diagonalization of the BC- hamiltonian Hbc which acts on 50(4) spinorial states \^e) 
like Hbc \^e) = E\^e) w hh Hbc = (dp — an< ^ w here the four matrices 7 M 

are the usual Dirac 4x4 matrices in the 4D euclidian space. By choosing the constant 



x By QHE, we mean a quantized conductivity of the system following from the existence of a discrete 
energy spectrum and a discrete filling factor due to the background field. 
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background fields like, 



•F^w — B £ nv3<i + ££i2nv , det T e ™ — B 2 x £ 2 

(1.1) 

2 

7=1 SU3 roots a 

breaking down the 50(4) symmetry of the euclidian M 4 down to SO (2) x SO (2), one 
ends with a diagonal form of the squared operator Hg C ) as well as remarkable factorized 
relations that allow to perform the explicit computation of the exact hamiltonian spectrum. 
Notice that, viewed from the 4D space time with SO (4) symmetry, the constants B and 
£, appearing in the above relations, are "magnetic" and "electric" components of T^v m 
respectively normal to the x x -x 2 and x 3 -x 4 real planes of R 4 . However, from a (1+4)- 
dimensional space time 2 with SO (1,4) isotropy symmetry, both of B and £ behave as 
magnetic components that couple to left and right handed fermions respectively; 
and so lead to QHE phenomenon in four dimensions. This 5D interpretation will be 
discussed with some details in section 3 and in conclusion. With the above diagonal 
choice of Tf$ , which contains as particular cases the singular limit det J-^J) = describing 
chiral configurations with solely or we find, amongst others and besides the 5D 
interpretation of the BC model, the two following features: 



(1) the energy spectrum the BC fermions in the constant background fields (1.1) is discrete 
provided the T e ™ tensor is non degenerate; that is as far as the product B x £ ^ 0. In this 
case, the energy spectrum is given by 



E^ m (ro, zu'^j = ±h\/ nw 2 + mm' 2 , n,m>0, (1-2) 

with the oscillator frequencies w 2 = , w' 2 = ^£i; and where Q = qe stand for the 
electric fractional charge of the quark u a (x) or the quark d a (x) living at each of the two 
valleys of the BC model. Clearly the appearance of the frequency w' is due to euclidian 
nature of the 4D lattice and to the underlying 5D interpretation where £ is thought of 
as an external "magnetic" field that couple to right handed *$>r. In the singular limit 
det = 0, for instance in case where £ — > 0, the above energy spectrum gets modified 
as 

I U2 i U2 

E±(m,k z ) = ±hmJn + ^-f^, 

V U7 

with hki being the momenta of ^/r along the i-th direction. These energies define a family 
of pairs of opposite paraboloids separated by the gaps AE n = 2-\J nw 2 + k 2 + k 2 , and touch 
on the fundamental state n = 0, k z = 0, k T = leading to a zero gap system; see also fig [j] 
for illustration. 



2 In 1+4 dimensions, the antisymmetric tensor Tmn has 10 components; 6 of them given by .F M „ are 
magnetic type and the other 4 ones given by are of electric type. 
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(2) The BC fermion's filling factor vbc of the anomalous quantum Hall effect induced by 
the background fields in the band energy < E+ m < E^ M , with N, M positive 

integers, is given by, 

(2JV + 1) (2M + 1) . s 

vbc = ky x N c x g L x g R x ± ^ '- (1.3) 

where ky is the number of Dirac valleys, which takes the value ky = 2 for minimally doubled 
fermions including BC and KW ones; ky = 16 for naive fermions and ky = 4 for staggered 
fermions; N c the quark's color number which is equal to 3; and gi = 9r = 2 the number of 

spin polarizations of the left and right handed fermions of the SO (4) spinor. Notice that 

SU (3) 

in the case where, in addition to = Be + Se^v, we also have 7^ along 

the two directions of the Cartan subalgebra of the SUc (3) gauge symmetry, BC fermions 
coupled to the background fields are described by 6 oscillators [2 for each color; i.e: for 
each ^> C L and with c = 1,2,3) and the above filling factor generalizes as follows, 

1 6 

4c = kyxg L xg R x-H (2N, t + 1) . (1.4) 

i=i 

The presentation is as follows. In section 2, we review briefly the anomalous quantum Hall 
effect (AQHE) of graphene; this section aims also to revisit useful aspects of graphene using 
relativistic field theory in (1 + 2)D ; and also to describe our approach on a simple system. 
In section 3, we study the minimally doublet fermions in a constant background field by 
focusing on BC fermions. In particular we study the BC model of 4D lattice QCD; first 
as a lattice field theory on the hyperdiamond; and second as a (1 + 4)D extension of the 
graphene near the Dirac points. In section 4, we study the algebra of the gauge covariant 
derivatives and its highest weight representations. These representations are used in section 
5 to determine the spectrum of the BC-hamiltonian in presence of background fields; and 
determine the filling factor of the AQHE of the BC fermions. We also give the relation 
with the spectral flow hamiltonian considered in (34|, [35], |{| p| pq| . Last section is devoted 
to conclusion and comments. In the appendix A and appendix B, we develop further the 
link between our study and the index theorem. 

2. AQHE in 3D relativistic systems: case of 2D graphene 

In this section, we study the anomalous quantum Hall effect of a relativistic fermionic 
system described by the (1 + 2) Dirac equation in a constant magnetic background field. 
This concerns a particular QED3 model where the electromagnetic field strength F^ v takes 
a constant value with direction normal to the 2D space surface where live the sheet of 
graphene. We assume that the magnitude of the external field B is bounded as follows 14 



Tesla < \B\ < 20 Tesla so that the Zeeman coupling can be ignored [40]. 



2.1 Dirac equation in 3D 

First consider the (1 + 2) space time M 1 ' 2 with local coordinates = (t,x,y), embedded 
in the usual (1 + 3) dimensional space time R 1 ' 3 parameterized by X M = (x^,z); then 



- 4 - 



focus on the Dirac equation of a fermionic particle, described by the complex field doublet 
ip a = ((f), x) living in M 1 ' 2 , in the presence of an external constant magnetic field B taken 
along the z-direction of M 1 ' 3 , 



2 2 



LL<W = o. (2-1) 

fi=0 b=l 

In this system of equations, the hermitian 2x2 matrix V = ia^D^ with iD^ = i (d^ — i\A^) 
and (iD^) + = iD^, is the gauged Dirac operator in (1 + 2) spacetime dimensions, cr° = I 2 
the unity matrix and a' 1 = (a 1 , a 2 ) the usual 2x2 Pauli matrices obeying the 2D Clifford 
algebra crV- 7 + a^a l = Ib^ as well as the SU (2) commutation relations type [a 1 , cr 2 ] = 2ia 3 
giving the ± chiralities of the component fields <p and %. Moreover, is the gauge potential 
which in the present case reads as \T illV x v where the constant field strength tensor 

T^ = -[D^D V ] (2.2) 

te 

is given by the magnetic field Beq^z along the z-direction. Here, the rank 4 tensor e a(U ,p 
is the completely antisymmetric SO (1,3) invariant tensor; and the restricted £0^3 = £\w 
is the antisymmetric SO (2) invariant 2x2 matrix. Substituting J 7 ^ = Be^, we see that 
the linear gauge potential lives in the x-y plane with components as given below, 



A 1 = +>§y, A = 0, 
A 2 = -fx, V^ = 0. 



(2.3) 



From the above relations, we learn also that [D\, D2] = i\B showing amongst others that B 
acts as a deformation parameter on the free spectrum of the Dirac operator; and the algebra 
of the gauge covariant derivatives, is up to a scale factor, the well known Heisenberg algebra 
of the harmonic oscillator; i.e [A, A^] = hi. Furthermore, factorizing the time dependence 
of the wave function by summing over all possible frequencies as 

/+00 jp 
|^e-^V|(z,y) (2.4) 

or equivalently by using discrete spectrum notation like 

+00 

r(t,x,y)= e- 1 n E - t ^(x,y), (2.5) 

n=— 00 

with E n = hui n ; and V'+| n | ( x -> u), ^-{nl ( x >2/) the wave functions associated with positive 
and negative frequencies; then using the fact that scalar potential A = 0, it is convenient 
to rewrite the Dirac equation of V'n as follows 



i(D 1 - iD 2 
i (D 1 + iD 2 ) 




(2.6) 



The solutions of this matrix equation give the possible frequencies oj n and the corresponding 
wave function doublets (4> n ,X n )i w ith the integer n refereing to the quantized values of 
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energy and momentum due to the background field B. To that end, acting on the above 
2x2 matrix equation ( |2,6| ) by the Dirac operator i ^ cr'-D/, we can bring it into the following 
diagonal one 

/ D_D+ 




\ D+D. 

where we have set D± = i (D\ ± 1D2) with the property (D±) + = D T . We also have 

^+=^( 2 i-f^) - U = X + l y, 1 = 1-9.-^9. ; 
D-=i(2l + §n) , u = x- iy , ^ = \l + hl y , K ' } 
together with the angular momentum L z around the z-axis reading like 

__ / d_ _ _d_\ _ 1 / d_ _ d_\ 
z \ du du) i \ dy ^ dx ) ' 

and acting on D± and D_D + as follows 

[L Z ,D±} = ±D± , [L z ,D-D + } = . (2.10) 

The relation [L z , D-D+] = shows that the wave functions carry, in addition to energy, 
also an angular momentum charge that we skip below for simplicity. Notice moreover that 
the D± operators satisfy the commutation relation [D-,D + ] = ^-B which, up on using 
the normalization 

D_ = A^B, D + = A^^B, (2.11) 

with positive B, can be put in the usual Heisenberg form 

[A,A]=I , [A,A)=[A^£]=0 , (2.12) 

describing the algebra of the quantum harmonic oscillator with A^ the creation operator, 
A the annihilation operator; and HI (with h set to one) the diagonal one that captures 
the energy quantum number n. The highest weight representations of this algebra are 
well known; they are given by semi infinite dimensional vector spaces generated by the 
normalized basis vectors {^± n ( x > u) = l£±n > } n >o based on the highest weight vector |£ > 
obeying the conditions 

A|£ >=0 , hl\t >=h\£ > , (2.13) 

solved generally like £ (u, u) ~ u l e~^ uu with I the value of the angular momentum which 
we set below to zero. The other states of the representation are given by successive appli- 
cation of the creation operator A^ as given below, 

At|£ n > =V^+l|C„+i> , n>0 

A\Z n > = V»K»-i> , TV = A^A 1 • ' 

with >= as well as 

M\Z n >= n \U >, <Zn\£m>= S nm- (2-15) 
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Substituting the operators D± in terms of their expression using A and A\ we can split 
the matrix equation ( |2.7| ) in two eigenstate equations as follows 



{A + A + 1) 4> n - 2if„ , (2 10) 

A+ A v - ^ v ' 

^ 71 A — 2e/3 An 

The wave functions of these relations are naturally solved by the basis states {|^ n >} by 
taking <p n and \n like 

(\Z n -\>* 



(2.17) 



describing two neighboring polarized particles. The energies follow from the relation be- 

9 E 2 

tween the squared frequencies io n = and the integer n. We have 



UJ 



2 



, E± = ±J^\n\ , (2.18) 



with E+ and E n respectively associated with the waves { x iV) an d ^-|n| ( X 'V) °f 



eq(2.5). Notice that for n = 0, we have a zero mode with negative chirality 

(2.19) 




2.2 Filling factor 

As was discovered in 1980, the Hall conductivity a xy of a 2D electron gas in a strong 

2 

transverse magnetic field develops plateaus at values quantized in units of i.e <7, 



2 . £ 2 

x v = V T 



with v = jjt- standing for the filling factor 44|. In the case of 2D graphene, the 



filling factor v gra ph of the QHE is given by the ratio of the number Np of particle states 
with respect the number of quantum fluxes = j s2 J~2 where the 2-form T% is given by 
Bdx A dy. We have, 

N f 

Vgraph ~ ~Tt~- (2.20) 

A direct way to get this factor is to compute the number Np per a unit quantum flux 
Nfj, = 1. For that purpose we proceed as follows: First recall that the ip + \ n \ and ip-\ n \ 
components of the expansion of the quantum wave function ip = ^ n ip n with ip n (t, x, y) = 
e~ tU)nt ip n (x, y) as in eq([E^) describe respectively particles and holes that live on the energy 
levels E+ and E~ of eq d2.18[ ). Second, the total number of Np polarized particles and Np 
holes, per unit quantum flux, that fill the band energy 

0<E%< E 2 N , E N <E n < E+, (2.21) 



with Ejj = dby \N\ is, because of symmetry, precisely 2Nf; and is given by the relation 

2N F = dtd 2 x\^(t,x,y)\ 2 (2.22) 
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together with the expansion 

+N 

iP(t,x,y) = J2e~ i " nt ^ n (x,y)- (2-23) 
—N 

Putting back ( |2.23| ) into ( |2.22j ) and using the orthogonality property J d 2 x ^J n (x, y) ip m (x, y) 
= <>n,m, the relation ( |2.22j ) reduces to 

+N 

2N F = 1 = (2JV + 1) . (2.24) 

-N 

From the above relations, we learn that the positive energy values are associated with 
the wave functions indexed by positive integers n > namely V>+ n x i u) = e ~ ^ Ent tp+ n ( x > v)'i 
they correspond, in the language of condensed matter physics, to the conducting band 
where live particle states. The negative energy values E~ are associated with the wave 
functions indexed by negative integers n = — \n\ < 0, V^m (t,x,y) = e~n Ent ip\ n \ (x,y); 
they are associated with the valence band where live the holes. In the particular case n = 0, 
both energies E^ vanish and the conducting and valence bands touch. So the fundamental 
state ipQ with zero energy would be filled by \ particle and \ hole; this property can be 
viewed as another statement of the chiral anomaly. 

Therefore the number of polarized particles filling the band energy < E^ < ^ 2!i ^ B \N\ is 
equal to ^V"+o- So, by taking into account the fact that electrons has two spin polarizations, 
the general expression of the filling factor v gra ph reads as follow: 

v graph = 2k v (n + ^ = AN + 2, (2.25) 

where ky stands for the number of Dirac valleys which is equal to 2 in the case of graphene. 

3. Borigi-Creutz fermions in background fields 

In this section, we extend known results on 2D graphene to BC fermions of J^D lattice QCD; 
but the analysis applies straightforwardly to KW and naive fermions. First, we recall the 
BC fermions and its two Dirac valleys as being a simple model for lattice QCD simulations. 
Then we study some remarkable properties of the Dirac equation in the background field 
Ql . 1| ) . The explicit computation of the spectrum of BC fermions in background fields will 
be given in next sections. 

3.1 BC model and 4D graphene 

Borigi-Creutz fermions is a simple four dimensional lattice QCD model for studying 
and simulating the interacting dynamics of the two light quarks up u a (x) and down d a (x) ; 
see ref.|7| for other specific features of BC fermions; in particular the explicit derivation of 
the SU (2) flavor 3 symmetry (u a , d a ). The BC model can be derived as a particular model 



3 To be precise the species in lattice fermions should be treated as point splitting fields as proposed in 
and used extensively in |l0[. 
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resulting from the following general euclidian lattice action 



Xi eL \/=i / 



(3.1) 



where L is the J^D hyperdiamond with lattice parameter a; \1/ Xi is a 4D euclidian Dirac 
spinor at the site x%\ 7^ the usual 4x4 Dirac matrices to be specified later; are five 
relative vectors parameterizing the first nearest neighbors; and where fi^ is a 4 x 5 tensor 
that reads like @, 



l+2i 


1 


1 


1 


\ 


4 


4 


1 


4 


1 


l+2i 


1 


1 




■1 


4 


4 


4 




1 


1 


l+2i 


1 




4 


4 


4 


4 




1 


1 


1 


1+2 




4 


1 


1 


4 



(3.2) 



Expanding the Dirac spinors in Fourier modes, we get Sbc = f ,t u ^w^BC^k where the 

J [2nj 

4x4 matrix T>bc is the 5C operator in the reciprocal space; it reads as follows, 
444 

V BC ~ iE^ Silla ^~ aE^ COSa ^ + aE rC ° Sa ^~ a 17 > (3-3) 



/*=1 



#*=1 



with 



r = i( 7 1 + 7 2 + 7 3 + 7 4 ) 



(3.4) 



From these expressions, one can check that the operator T>bc has two zer ° modes given 
by the two following four component wave vectors k^, 



(i) : K ={kiMMM) =(0,0,0,0) 
(m) : K> =(h,k 2 ,k 3 ,h) = • 

The first zero mode is associated with the quark u a (k); and the second with the quark 
d a (k) , these fields are just the Fourier transform of u a (x) and d a (x) respectively. Ex- 



(3.5) 



panding eq (|3.3|) near the zero modes (|3.5|-i) and ( |3.5| -ii), one gets at the first order in 

A A 

d 



V B c - E 7 m ^m ~* E 



r 



(3.6) 



which is precisely the Dirac operator of a free fermion. As such, near each one of the two 
Dirac points the BC model reads in the continuous space as follows 



S BC = I d 4 x^+(x) 7 < 
Jr 4 



(3.7) 



with 'A? (x) standing for the quarks it (x) and d (x) respectively given by integrals over 
momentum hyper-balls B centered at K and K' as reported below; see for technical 
details, 

^W = E/B K ^ ik(x ^(y) , 

d ^) = tk Kl ^^-y)d{y) . (3 - 8) 
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Notice that the way this the above relation is a naive definition of the flavors; a more 
concise way to do is that done in ||]. Although not used here below, this method is based 
on point splitting fields in the reciprocal space. For later use, it is interesting to notice as 
well the following remarkable properties: 

First, one has to think about the field ^ (x) as the zero mode of the hamiltonian equation 



{x\H bc \^e) = E^ E (x) , {x\H B c\Vc 







with hamiltonian 4 



H B c = 7 A 



d 



(3.9) 



(3.10) 



To interpret eq([3.9D as an equation of motion, one has to modify the action Sbc ( |3.9| ) by 
adding an extra term of the form M- (x) ^ e ( x )- 

Second, since the parameter E can take any real value, we will also think about the Dirac 
spinor ^e (%) as describing modes of a spinor field living on a 5D space time as given below 



*(x,t) 



dE 
2ttE 



e* Et * E (x), 



(3.11) 



so that the J^D action looks as resulting from the following five dimensional one 



S5D 



j 



dtd A * 



(3.12) 



where t stands for the five dimensional coordinate; and which could be interpreted as a 
"time" coordinate encoding non zero energy deformations near the Dirac points. This 5D 
interpretation will be developed with further details later on; it will be also linked with the 
so called spectral flow hamiltonian given in sub-section 5.3 and appendix B. 

3.2 BC fermions in background fields 

We start by giving some useful tools on spinors in the euclidian space time M 4 ; then we 
describe the Dirac equation of the BC fermions in a generic background field J-^. 

3.2.1 From SO (1,3) to SO (4) and then to 50(1,4) 

First, recall that the quarks are confined relativistic particles described by four dimensional 
space time fields u a (x) and d a (x) transforming as complex 4-component Dirac spinors, 



*a (x) 



/*i(x)\ 

* 2 (x; 
* 3 (x; 

\*4 (X)/ 



(3.13) 



4 Notice that Hbc has been defined as an invariant under the SO (4) symmetry that follows from the 
continuum limit of the hyperdiamond. If one insists on defining the hamiltonian in terms the translation 
generator along the 4-th direction, one has to break SO (4) down to SO (3); and ends amongst others either 
with detTfiv — or with a quantization condition on the background fields B and £. 



- 10 - 



These particles carry several charges, in particular the spin charge of the SO (1,3) space 
time symmetry as well as the electric and color charges respectively given by the U em (1) 
and the SU C (3) gauge symmetries. The color charge is understood in eq( 3.13|) ; since each 
component \I> a should be read as ^ c a with c = 1, 2, 3 the color index; i.e: 



* a 



(3.14) 



Recall also that in 4D lattice QCD, one uses euclidian time rather than the Minkowski 
one, so the standard SO (1,3) Lorentz symmetry gets mapped to SO (4) which is locally 
isometric to the product 

SU L (2) x SU R (2) . (3.15) 
k.r — k^x 4 get replaced by the SO (4) ones namely 



So SO (1,3) scalars type k^x^ 
k.r + k^x 4 ; similarly the SO (1,3) semi-norms k 
k 2 + k\ = C 2 leading to k 4 = ±\/~C? 



k\ 



£ 2 transform into positive norms 
k 2 with imaginary values in the case where £ = 0. 
Moreover, the four component SO (4) Dirac spinor ty a ( |3.13| ) is essentially made by the 
direct sum of the two SU (2) spinors 



R, 



(3.16) 



where undotted fermions refer to SUl (2) and dotted ones to SUr (2). Each one of these 
Weyl spinors have two components as described in the following table, 



symmetry groups 


SU L (2) x SU R (2) 


~ SO (4) 


S!7 C (3) 


representations 


fro) 




(i,o)e(o,±) 


3 


fields 




L 


•-(£) 







(3.17) 



The fields (x) and | A (x) are just the left handed and right handed quarks ^ l (x) and 
$ (x) ; they carry opposite chiralities under 75 and live on the four dimensional Euclidean 
space time M 4 parameterized by the local coordinates x = (x, y, z, r) with diagonal metric 
b^ v = diag (+ + ++). These fermions have integer color charges; but fractional electric 
ones given by Q u = y and = — §. 

Moreover, in lattice QCD the BC- fermions ^ (k) living at the Dirac points ( |3,5] ) have 
zero energy as they are zero modes of the Dirac operator. Therefore non zero energy ex- 
citations induced by deformations near the two Dirac points may be interpreted in terms 
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of the 5-dimensional direction t introduced above; so that the space with local coordinates 
x = (x,y,z,r) gets promoted to a (1+4) space time with coordinates 



(x;t) = (x,y,z,r;t) . 

As such the four component field ^ (x) gets also promoted to & (x;t) = Sfr (x,y, z,T]t) 
whose expansion with respect to the fifth coordinate t is given by eq( |3.11| ). 

3.2.2 Dirac equation coupled to 

The Dirac equation describing the dynamics of the four component Dirac fermion ^ c a (x) 
coupled to the background vector gauge potential (x) is given by the following system 
of four coupled equations 



V) cd 



(3.18) 



where (^4^) cd = Yli Af/Tj d with the generators of G = SU C (3) x U em (1) and go the gauge 
coupling constant. Below, we will restrict to the particular case where A^ em{1) / and 



SU C (3) 

!/=m(l) 



0; the case with T[ 



■SU C (3) 



/ and A 



SU C (3) 



/IV 



7^ will be discussed in conclusion section. In the case 



0, the above equation reduces to, 



da - i—A^ 

c 



(3.19) 



where now A^ stands for the electromagnetic gauge potential A^ em . Notice also the two 
following things: 

First, by help of the constraint equation = expressing the fact that is a constant, 
it is not difficult to check that the electromagnetic potential vector A^ is linear in space 
time coordinate positions as follows, 



A, 



2 -> liV'L i 



d»A L 



"-t = W U = \Tr (7) = 0, 



(3.20) 



where we have dropped out the irrelevant integration constants. Generally speaking, the 
antisymmetric tensor = d^A^ — d^A^ capture six real degrees of freedom, three com- 
ponents for the "magnetic" field B and three components of the "electric" field £, 



/a' 



( +B 3 -B 2 +£i \ 
-B 3 +Bi +E 2 
+B 2 -B x +E 3 

V -£ x -E 2 -£ 3 / 



det T, 



(e.bY 



(3.21) 



From the 5-dimensional view, this matrix should be reads as a 4x4 submatrix of the 
following 5x5 antisymmetric one, 



MN 




(3.22) 
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where now the J-^'s are the 4 components of the electric field in 5-D; and Ta V the 6 com- 
ponents of the magnetic tensor. Below, we will focuss on particular situations where J-^ v 



has two degrees of freedom as in eq(l.l). This corresponds to choose the four components 
of the electromagnetic gauge potential as follows: 



Ai 



2 »> 



At 



2 *"> 



(3.23) 



leading to the following constant electromagnetic field strength, 

/ —B \ 



+B 
-5 
\ +£ o y 



(3.24) 



with det J-^j, = ,B 2 x £ 2 . In this relation, B and £ are a priori arbitrary real numbers; but 
to have a quantum Hall effect, they have to be large enough so that interacting energy is 
much greater with respect to kinetic energy. 

The second comment concerns the four 7^'s of ( |3, 19 ); these are the J^D euclidian 4x4 
Dirac matrices obeying the usual J^D Clifford algebra, 



(3.25) 



with I 4 the 4 x 4 identity matrix. Notice that the SO (4) ~ SU L (2) x SU R (2) has 6 
generators, 3 of them generate SUl{2) and the three others generate SUr(2). In the 
spinor representation, these are given by the Pauli matrices t 1 and a % with generic 2x2 
matrix representation is as follows, 



(3.26) 



with Q l standing for both t % and a % . The matrices g , g 2 satisfy the 2D Clifford algebra 
g % gi + gig 1 = 25 13 12 and the SU (2) symmetry bracket [^g 2 ] = 2ig ?J . Notice moreover 
that the euclidian 7^ matrices can be realized in terms of tensor products of r l and a % as 
follows, 






7 4 = r 1 ® a 4 



7 5 = r 3 ® a 4 



T° = r 4 ® a 4 



(3.27) 



with t 4 , a 4 = I2 = I and 7 s = 7 X 7 2 7 3 7 4 . More explicitly, we have 



7 



1(7 



-ia h 

k 







7 



I 
/ 




^0 




(3.28) 



The anticommutator of the 7^ matrices verify the J^D Clifford algebra ( p. 25 ) and their 
commutators 7^] give precisely the 6 generators of the spinorial representation of the 
SO (4) symmetry, 



7 4 7 i 



7 J 7 * 

y 7 4 



2ie^ k (1 



2ie 



123 



)(T*) 



(3.29) 
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where £^ k is the completely antisymmetric 3D Levi-Civita tensor. 

Now, we turn to determine the spectrum of BC fermions in background fields by solving 
eq( 3.19| ); but before that it is interesting to study the properties of the non commutative 
algebra of the covariant derivatives; then turn back to compute the spectrum. 



4. The Lie algebra of the gauge covariant derivatives 

As we will see, the algebra of the D„ covariant derivatives in constant background fields 
has much to do with Heisenberg algebra of quantum harmonic oscillators. Because of the 
tensor structure of J 7 ^, we distinguish two basic cases: first the case of two uncoupled 
quantum oscillators associated with the choice Q3.24 ) and detJ 7 ^ 7^ 0. Then we con- 
sider the general case ( 3.21| ) with det J 7 ^ = (£.£>) 7^ describing two coupled quantum 
harmonic oscillators. 



4.1 uncoupled case: = Be ^4 + £e\ 2fll/ 

Generally, the four gauge covariant derivatives D\, D 2 , D 3 , D4 satisfy the generic commu- 
tation relations 

Q 



-i—T, 



but because of the choice ( |3.24 ) of the background fields, they reduce to, 

[D 1 ,D 2 ] = i9£ , [D 3 ,Di] = i&- , 



(4.1) 



(4.2) 



and all others vanishing. These are very special relations; first because their right hand 
sides are non zero constant numbers showing amongst others that the large limits of B and 
£ induce a non commutative geometry in the euclidian space, 



with Q^ v given by 



Q 



( -£ 3 £ 2 -BA 



(4.3) 



->av 



1 



£.B 



(4.4) 



£3 -£ x -B 2 
-£ 2 £ x -B 3 
\B 1 B 2 B 3 / 

Second, the four covariant derivatives organize in 2 + 2 capturing a complex structure as 
follows, 



Di ~ iD 2 = I + §f u , D 1+ iD 2 = ^-^ 
D 3 -iD, = §, + ^v , D 3 + iD± = ™-%v , 



2d _ QB 



II 



(4.5) 



together with similar relations for their adjoint conjugates. In getting eq fl4.5| ), we have used 
the explicit expressions D\ = d\ — i^-y, D 2 = d 2 + i^r-x together with similar relations 
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for D 3 , D 4 ; and we have set 



u = x + iy 



A. - I ( A. _ j£. 

' du ~ 2 \ dx L dy 



v = z + ir , & = 



(4.6) 



To study the representations of eqs( |4.2| ), one has to specify the sign of Q— and As- 
suming ^— > 0, > for simplicity and which means that Q, B, £ have the same sign; 
and setting 



i(D 1 -iD 2 )=AJ®£ , i(D 1+ iD 2 ) = A\ 2 -^ 



i(D 3 -iD 4 ) = C^ 2 Sf , i(D 3 + iD 4 ) = C^ 2 ^ , 
the commutation relations (14.21) read also as 



(4.7) 



[A,Ai]=I , [C,Ct]=/ 
[A,C] = , [A,&]=0 



(4.8) 



These relations (4.8) show that the Dirac fermion in the background field fl3.24| ) de- 
scribe a priori two quantum harmonic oscillators with oscillation frequencies 



w 



2QB 



— 



2Q£ 



(4.9) 



The operators A^A and C^C giving the number of energy excitations in w and w' units 
respectively read in terms of the gauge covariant derivatives as follows 



225 AU = (D 1 ) 2 + (D 2 ) 2 + i[D 1 ,D 2 ] 
= (Dr) 2 + (D 2 ) 2 - Of 



(4.10) 



(D 3 f + (D 4 ) 2 + i[D 3 ,D 4 ) 



(4.11) 



and similarly 

= (D 3 ) 2 + (D 4 ) 2 -<f 

Using the space time variables (x,y,z,r) and their derivatives d^, these operators read 
more explicitly like 



(£l) 2 + Wj 



a 2 



dx 2 dy 
id 2 



+ 



dudu 



+ 



QB 

2c 



QB 

2c 



UU + 



(x 2 + y 2 ) + 



QB 
2c 



C'xy > 



QB 

2c 



r - 



(4.12) 



and 



(D 3 ) 2 + (D 4 f 



a 2 _ a 2 

dz 2 dr 2 



+ 



QS_ 

2c 



QL 

2c 



r 

'-'ZT 1 



id 2 



+ 



QS_ 

2c 



(vv) + 



Q£_ 

2c 



r - 

>~-vv > 



(4.13) 
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where one recognizes the harmonic potential energy induced by the electromagnetic flux 
and the angular momentum operator components C xy and C zr 



describing respectively rotations around of the oscillations of the particles the x-y and z-r 
planes of the 4D space time. Notice that the hermitian operators C xy and C ZT are nothing 
but the charge operators of the abelian Ul (1) and Ur(1) sub-symmetries of the SUl (2) 



and SUr (2) respectively. Indeed, using the complex variables (4.6) and putting back into 



( |414|) , we find, 



acting on the monomials u ni u n2 v mi v m2 as follows 

[£ xy ,u ni u n2 v mi v m2 } = (n 1 -n 2 )u ni u n2 v mi v m2 

(4.16) 

[C ZT , u™ 1 u n2 v mi v" 12 } = (m 1 -m 2 )u ni u n2 v mi v m2 . 



Notice that for the particular cases n\ = n 2 and mi = m^, there is no Ul (1) nor Un (1) 
charges. So wave functions \I/ a = ^f a (u, v, u, v) that depend only on uu and vv have no 
angular momenta. 

Moreover, using the expression of the matrices 7^, we can write this matrix operator 
{HBc) a j3 = ^ap^n m terms of the creations operators , and the annihilation ones 
A, C as follows: 



(H BC ) 



( 


w'C wA 
\wA^ -w'C* 



-w'C* -wA\ 
-wA* w'C 



/ 



(4.17) 



with H BC = Hbc- Moreover, using the commutation properties CA^ = A^C and CA 



AC and their adjoint conjugates following from the choice ( p. 24 ), the squared BC hamil- 
tonian H BC reads as follows, 



(w 2 AA^ + w' 2 C^C 

w 2 A^A + w l2 CC^ 

w 2 AA^ + w l2 CC^ 









\ 



V 















w 2 A^A + w' 2 C^C j 



which allow to determine spectrum of the Borigi-Creutz spectrum; this will be done after 
studying the coupled case. 
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4.2 coupled case: = B^e^j + £ M 



In the case where the 6 components and £[ 4l l, = 1,2,3) of the background field 



T^v are non zero as in ( |3,21| ), the previous two Heisenberg algebras (4.8) get coupled. 
Indeed, using the relation [L> M , D v \ = —i^-J 7 ^ and eq(|4.7|), we find in addition to 



[A At] 



[c,ct] 



the following couplings 



[AC] =A, [A,Ct] 
[Act] =Cj [a,c] 



-A, 

-c, 



where we have set, 



A 



l 



2VBS 
1 



[(^14 + ^23) + * (^13 - ^24)] 
[(^23 - -Fu) + * (^13 + -^24)] 



(4.18) 



(4.19) 



(4.20) 



Using the BC hamiltonian ( 4.17 ), we find that its square H% c is given by, 



(w 2 AA^ + w' 2 C^C -zjw'C \ 

-W( ro 2 AtA + ro /2 CC7t 

w 2 AA^ + w' 2 CC^ -ww'\ 

\ - wa ro 2 A.4 + ro /2 ctc y 



and is no longer diagonal. Moreover, the operator numbers A' A and C^C do no longer 
commute; a difficulty that may be overcome by treading A and £ as small perturbation 
parameters. 



5. More on AQHE in BC model 

First we compute the spectrum of BC fermions in the constant background fields; then 
we determine the filling factor vbc of the associated quantum Hall effect. After that, we 



study the link with spectral flow method of [34, 



5.1 Spectrum of BC Hamiltonian Hbc 

Starting from the BC fermion's equation Hbc = E \ ^e) , then substituting the matrices 



7^ by their expressions (|3.28j ) and the Dirac fermion ty a by its two Weyl spinors (*0 a ,£<j), 
we get 



O 
0t 



E 



(5.1) 
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where the O and 0t operators are given by 

n ■(-^'C ] -wA\ f (w'C wA \ 

These operators can be also written in the following condensed form 

3 3 

O = a 4 D 4 + i ^ ° k D, ] = a A D A - i^^D, (5.3) 

k=l k=l 

which let understand that solutions may be as well obtained by using quaternions. To deal 
with eqs( |5.lD , we will use an explicit method relying on acting on its both sides of by Hbc- 
We get the following set of equations 

ooH = E 2 ^ , 

O^Ol =E 2 l , 1 j 

showing that the two Weyl spinors and \ h can be treated separately; thank to the choice 
Q3.24 ). This property was also expected as there is no mass term that couples the left and 



right components of the Dirac fermion. To solve ( |5.4| ), we use the following expression of 
00t and 0*0, 

UU ~ [ o w' 2 [CtC + 1] + zu 2 A^A J ' (5 ' 5j 



and 



n \n - ( w ' 2 i° 1c + !] + w2 i AlA + 1 \ 

w' 2 C^C + w 2 AU 



(5.6) 



Notice that these matrix operators don't commute; they obey the commutation relation 

OO^-O^O = -w' 2 ^ . (5.7) 

Then putting the expression of OOt and 0t0 back into eq( |5.4|) , we obtain the following 
system of 4 equations which allow to determine the 4 component of the Dirac spinor 

[w' 2 C^C + w 2 [A^A + 1] ] = E 2 ip 1 , 

[w 12 [C7tc + l] + W 2 A^A\ ip 2 = E 2 ip 2 , 
[ w ' 2 [C^C + l]+w 2 [A^A+l]\l l =E% , (5 ' 8) 



[ro ,2 CtC + w 2 A^A] l 2 = E 2 l 



2 
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To solve these relations, it is useful to introduce the wave function basis {|i? n ,m)} of two 
uncoupled quantum harmonic oscillators. These wave functions are factorized as 

A n , m (u, u, v, v) = n (u, u) x 6' rn (v, v) , n,m>0, (5.9) 

with the constraint relations 

Ax6 (u,u) = (f + 6 (u,u) = 



cxo> (v,v) = (% + q§v) e (v,v) = o 

solved as follows 



(5.10) 



e (u,u)=M (u)e-^ uU , tf Q {v,v)=N( i {v)e-%™ , (5.11) 

where the complex functions Wo (u) and Mq (v) obey the holomorphic conditions 9J ^q^ = 
and d ^Q^ = 0; and are generally fixed by the charge of the angular momentum operators 
and the normalization condition of the wave functions. The others waves functions 9 n (u, u) 
and 0' m (v, v) are obtained as usual by applying the creation operators. 

e n (u,u) ~ ^-^ u y eo ( u ,u) , 



Because of the algebraic properties [A^ A, (^)™] = n (A^ n and [C^C, (C^)' 
the |A njjn ) states obey amongst others the following properties 



m 



A^A | A. 
CtC | A 
A^AC^C\A n 



n,m/ 
rn ) 



^ |A ni?Tt ) 
nm |A„ )m ) 



Taking the two Weyl spinors ip a and £ d like 



with remarkable chiral zero mode 



An i 772—1 

An-1, m-l 

\ A Hj m 



(ct) m , 



(5.13) 



(5.14) 



(5.15) 



/ \ 



V A o,o J 

and putting back into OO^tp a = ^,m^« an d Ot0£ d = E 2 m ^ h , we obtain the following 
the energy spectrum, 

E± m (w,w') = ±hy / nm 2 +mw' 2 . (5.16) 
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This relation is comparable to fl2.18| ) of graphene namely £ ^ D -9 ra P hene _ -)_y'2^B jy ^ ce 
also that in the case of the background fields B ^ 0; but £ = 0, the tensor is degenerate 
since det J-^ v and the algebras of the covariant derivatives (4.2) reduces to 

OB d d 

[ Dl ,D 2 ]=i^, D 3 = -, D 4 = -, (5.17) 

and all remanning others which are just commutation relations vanish. Therefore one has 
only one quantum harmonic oscillators with frequency w whose energy spectrum is given 

by, 

Et (w, k) = ±h\Jnw 2 + k 2 (5.18) 

with k 2 = fc| + k 2 and p 3 = hk 3 , p± = Hk^ are the momenta along the z- and r- directions. 
Notice that for the case Pz = Pa = 0, the energy levels are given by (w, k = 0) = 
±hm^/n with n a positive integer. Notice also that for the fundamental state n = 0, the 
energies vanish Eq (vj, k = 0) = Eq (w, k = 0) = 0; and the valence and conducting bands 
touch; for an illustration see fig |]. 



n=2 n=1 




Figure 1: the first levels of the energy spectrum {ru,k) = ±hvj\Jn + J^- of BC fcrmions in 
the background field T^ IV = Be with det T^ v = 0. For n = 0, the conducting and valence bands 
touch at the Dirac point k = 0. 



5.2 Filling factor vbc of BC fermions in background fields 

Before computing vbc, we start by recalling that the wave function ^ (x, t) of the BC 
fermions in the background fields ( |3.24| ) is given by 

00 

*(x,t)= e"^"' mi *n, m (x) (5.19) 

n,m=— 00 

with wave modes ^> n ,m = (^n,m> £n,m) an d energies E n ^ m respectively given by eqs( |5.14| ) 
and ( |5.16| ). Since the waves $„ im form a basis and using probability interpretation of 
\I/ (x, t) we also have, 

/ d 4 xV* n m (x) # <m , (x) = 5 n ^5 m ^ m > . (5.20) 

JR 4 
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These relations can be also derived by help of well known orthogonality of the waves 
functions 9 n (x, y) and m (z, r) of the quantum harmonic oscillators. Now, using these 
relations, one can compute the filling factor vbc of the AQHE of the Borgi-Creutz fermions 
in the constant background fields B and £. This factor is given as usual by j^; that is the 
ratio of the number Mf of particle states with respect the number of quantum fluxes Mq. 
The latter number reads as follows: 



^ = WZ I dx " A dy v T^. (5-21) 
Is 2 



2tt 



and, by using eq(3.24), can be split as Mb = f S 2 Bdx A dy and Ms = f S 2 £dz A dr. For 



a unit flux M^ = 1, the number Np of polarized particles that fill the band energy < 

n,m — N,M 



< E 2 N M is given by the relation 



M F = - I dtd 4 x\V(x,t)\ z (5.22) 



Substituting * (x, t) by eq( ^.19|) and using the orthogonality properties, we get 

AV = (2iV+1) 2 < 2M + 1) . (5.23) 

In the particular case N = M = 0, the energies Eq and Eq vanish; and the conducting 
and valence coincide; so the fundamental state ^o,o with zero energy is filled by ^ particle 
and ^ hole. The number of polarized particles filling the band energy < E+ m < E~^ M 
is equal to Mf- Therefore, taking into account the fact that a BC fermion has four spin 
polarizations; and each fermion has three colors, the general expression of the filling factor 
vbc reads as follow: 

, (IN + 1) (2M + 1) 
v BC = 12k v ± (5.24) 

where ky stands for the number of Dirac valleys which, for BC fermions, is equal to 2 as 
in 2D graphene. 

5.3 Link with spectral flow and topological index 

In this subsection, we want to comment on the relation between the hamiltonian, that we 
have used in this study and which has been motivated from 2D graphene, and the so called 



spectral flow hamiltonian considered recently in [34] in connection with the theoretical 
foundation for the index theorem for fermions on lattices. More analysis on this matter 
and other issues are further developed in the appendices A and B. 

In studying the AQHE in 4-dimensions, we used the equation fl3.19|) describing the dynam- 
ics of a fermion near a generic Dirac point corresponding to a zero mode solving = 0. 
This equation, that we put in the following simple form, 

(D-E)V = 0, (5.25) 

depends on the spectral parameter E, which in relativistic theory, has a dimension of mass; 
it captures deformations away from the Dirac point and allows to open the gap between 
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Figure 2: On left the Dirac point where conducting and valence bands touch. On right the gap 
between conducting and valence bands induced by the mass term. 



the conducting and the valence bands as shown on fig y. Recall that in eq(5.25), the Dirac 
operator in the background field J-^ u as in eq( |3.24| ) is an antihermitian operator given by 

D = -fD^ = 7^ (fy - iA„) . (5.26) 

Notice that the corresponding hamiltonian used in this study namely H (E) = {D — E) 
looks quite similar to the one used in the method of spectral flow of [34]. There the spectral 
flow hamiltonian is given by the following 

H sp (m) = 7 5 (D - m) (5.27) 

with 

7 5 (£>-m)1' = A(m)¥, (5.28) 

and the eigenvalues A (m) depending on the spectral parameter m. Clearly the hamiltonian 
H (E) and the spectral one H S p (rn) are close cousins and are related as follows, 

H sp (m) = 7 5 # M . (5.29) 

The trick of multiplication by 7 5 has two nice and remarkable effects; it makes H sp (m) 
hermitian and allows to lift the degeneracy of energies with opposite chiralities. Indeed, 
by splitting the Dirac spinor ^ into its two chiral components *$>± like 

*+=UJ ' *~ = (x) ' ^± = ±*± (5-30) 

and restricting the study to the Dirac zero modes Dty = 0, the eigenvalue equation of the 
spectral flow hamiltonian reads as follows 

H sp (m) ¥± = =Fm^ ± , (5.31) 

and captures perfectly the data of fig [2] contrary to H (E) which fails in this matter since 
we have 

H (m)^>± = -m-S>± . (5.32) 
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It happens that the spectral flow hamiltonian H sp (m) is a powerful ingredient as it leads 
to exactly the topological index Ind(D) of the Dirac operator. Following [34, and the 
analysis developed in the appendices A and B, one can show that the spectrum of H sp (m) 
is, up to an irrelevant sign, the topological index 



Tr {^e- tD ) =N + -N_. 



(5.33) 



In this relation, N + is the number of zero modes with positive chirality; and iV_ is the 
number of zero modes with negative chirality. The difference N + — iV_ is equal to the 
topological charge Qtop of the gauge field; i.e: N + — iV_ = Qtop- For the case of the 
hamiltonian H(m), the corresponding quantity (restricted to the zero modes) is given by 



Tr (e~ tD ) | zero modes = N + + AL 



(5.34) 



it is not a topological quantity; it gets extra contribution from non zero modes. We end 
by noting that in the present 4D case, we have 



quarks q 




n q - 


degeneracy 


N q+ 


N q - 


u (x) 





1 


1 Qx 1 1 Qz I 





\Qx\ \Qz\ 


d{x) 





1 


1 Qx 1 1 Qz I 





\Qx\ \Q Z \ 



(5.35) 



where the integer \Q X \ \Q Z \ is the degeneracy of the zero modes; see eq( |7.36| - |7.37[) of ap- 
pendix A. The index theorem reads therefore N + — iV_ = (N u+ + Nj + ) — (N u - + N^-) 
= 2\Q X \ \Q Z \ in agreement with results of ||. More comments and technical details are 
reported in appendices A and B. 



6. Conclusion and comments 

In this paper, we have studied the anomalous quantum Hall effect in 4D lattice QCD with a 
special focus on Borigi-Creutz fermions in presence of a constant background field strength 
T^v Recall that BC fermions were proposed to simulate the dynamics of light quarks using 
four dimensional lattice QCD and where shown to be intimately linked the J^D extension 
of graphene on hyperdiamond. We have shown that in the neighborhood of the two Dirac 
valleys (|3.5|), the BC model has a 5D field theoretical interpretation with the following 
features: 



(1) the background fields given by antisymmetric tensor J 7 ^ appear as the 4x4 submatrix 
of a 5D field strength Fmn having 4 + 6 components distributed like 



MN 



•F iw -F u5 



In this 5D interpretation, the four components Tpx> describing a 5D electric field are equal 
to zero (7^5 = 0); while the remaining 6 components Tnv> which describe the components 
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of a 5D magnetic field, have been used to study the AQHE of BC fermions. This property 
explains also the behavior of the background fields B and £ as a magnetic fields respectively 
rotating left and right fermions and and leads to AQHE with cyclotron frequencies 



(2) the hamiltonian Hbc is precisely given by the euclidian J^D Dirac operator J2^=i i-l^D^ 
with = (dfj, — ii^TjxvX 1 ^ capturing interactions between BC fermions and the " 5D 
magnetic" tensor J-^ u . Because of the non zero flux induced by the background field J-^v, 
the commuting algebra of the four flat space translations turns into the non commutative 
[Dp,D v ] ~ iJ 7 ^, [D^Tvp] = 0. In the case where is chosen as Be + £ei2nv, the 
algebra of the D^s splits into two uncoupled Heisenberg ones, 

[Di,D 2 ] = , [D 3 ,D A ] = -iQ§I , 

whose spectrum as well as the associated AQHE have been explicitly studied in this paper. 
This splitting teaches us moreover that [D^,D U ] ~ iT^v and describes in general two 
coupled quantum harmonic oscillators which can be studied by using perturbation theory 
if the following conditions hold 

1^14 + ^2 3 | 2 + 1^13 " ^24| 2 « 2 1^12^341 , 
1^14 " ^2 3 | 2 + 1^13 + ^24| 2 « 2 1^12^341 • 

We conclude this work by giving a comment on the extension of the result given in this 
study to the case where the background gauge potential is valued in the Cartan subalgebra 
of SU C (3) x U em (l), 

U 2 (1) x U em (1) C SU C (3) x U em (1) . 

This subalgebra is a 3- dimensional abelian algebra generated, in addition to ^f/ em (i) = 
Q em I, by the two diagonal hermitian 3x3 matrices hi and Ji2 of SU C (3) which read in the 
Gell-Mann vector basis as follows 







(1 o\ 


;:: 


H 


0-10 


\oo 1 ) 




V / 




In presence of non zero constant background fields T^y 1 ^ 7^ 0, J-fy 7^ 0, J-ff 7^ 0, the 
gauge covariant derivatives reads as follows, 

D-d - x u (^F^ + i —hoF { - 2) + i -^F i - em) \ 

obeying the non commutative commutation relations [D M , D u ] ~ iT^ but now with T^y = 
^2 j hiT* that reads in matrix notation like 

/4t m) +4' ) +4' ) 

n jr( em ) _ + jM n 

jffl - j 
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Moreover, because the /ij's are diagonal matrices; each gauge covariant derivative splits 



into 3 components , Dp , Djp as given below, 



D (l) 


= fy- 


2c x 


K xern-> ufj, 


+ ^ + 3^ 


D m 


= fy- 


2c 


f n Jem) 
mem'' v/i 


- ^ + ^ 




= fy- 


2c x 


f n rH 


- 2*7$) 



and one is left with £ quantum harmonic oscillators whose uncoupled realization is given 
by the choice 

In this case the filling factor of the AQHE of BC fermions induced by these J-^s is given 



by eq(1.4). In the end, it would be interesting to study the analogue of the Zeeman effect of 
2D graphene in 4D lattice QCD by using models like the BC fermions we have considered 
in this study. 

We close this conclusion by two more comments. The first comment concerns the relation 
between the our results and the index theorem. This matter has been described in a 
condensed form throughout the paper; in particular in subsection 5.3. It is developed 
with details in appendices A and B. The second comment concerns QHE in lattice QCD as 
whole; the present study focused in BC fermions should be viewed as a first step which itself 
need further investigations; several issues like the explicit breaking of discrete symmetries 
in BC fermions flotl as well as spin Hall effect and the associated topological indices [50, 51 
haven't been addressed here. 



7. Appendix A: Index of Dirac operator 

In this section, we describe the index theorem of the Dirac operator D in background 
gauge configuration. This index provides the relationship between physical properties and 
the topology of the space in which live the fermions. First, we describe the index in 2- 
dimensions; both in continuum and lattice QFTs. Then, we consider the main lines of 
the index in 4-dimensions and its relationship with the so called spectral flow hamiltonian. 
After that, we make an heuristic interpretation of the filling factors obtained in this study 
in terms of the index of the Dirac operator. 

7.1 case of 2D graphene 

We begin by describing the topological index theorem of the Dirac operator in 2-dimensions. 
Then we apply the construction to the case of 2D graphene. 

7.1.1 the index theorem: continuous limit 

In the case of fermions living on a 2-dimensional surface in presence of external gauge 
fields, the index theorem gives a relation between zero modes of the gauged Dirac operator 
V = a^pDp, with gauge covariant derivative D^ = (d^ — iA^), and the flux Q = ^ J s Ti 
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going through S (J-2 = \j-^ v dx^ A dx u ). Writing this Dirac operator as 2 x 2 matrix T> as 
follows, 

V=( 0D A, D + = D 1 + ,D 2 (D+) , = _ _ 

yD_ J ' D_ = D X -iD 2 ' v +y v 7 

and the flux Q through 5 like Q = "ff" with I? = ^e 3 ^^^^, the index of the Dirac operator 
reads explicitly in terms of the number Nq {N^) of zero modes with positive (negative) 
chirality and the magnetic flux as 

N+ - AT- = Q. (7.2) 

This index is defined like Ind {T>) = Tr [a^e~ tXD ^ with t some spectral parameter and 
leads to the above equality which can be got by computing Ind (V) in two different ways 
and equate the two results. Let us describe rapidly how this is done. 
(1) use the property that the Dirac operator T> and its square T> 2 have the same number 
of zero modes. Since T> 2 has a diagonal form 

V 2 = r } ^V ii V u - % -[a^cj v \T iiV (7.3) 
as clearly seen on the matrix representation 

V 2 =( D+D - ° I , V 2 ^> = E 2 ^> , (7.4) 




one can do explicit calculations by using this property. Let show rapidly how this works. 
From above equation, we learn that the operators D+D- and D-D+ have the same non 
zero eigenvalues £ / 0; this means that they have the number of non zero modes and 

N+ = N~, n^O, (7.5) 

with ± refereing to positive and negative chiralities. Denoting by ^ = \P— ) the 2D 
spinor with and standing for the positive and negative chiralities; and assuming 
that 

D + D-t> + = \^ + , A/0 (7.6) 

then is an eigenstate of D_D + with the same eigenvalue; but with opposite chirality. 

This feature follows directly by multiplying both sides of eq (|7.6| ) by D_ . However in the 
case where A = 0; that is D + D_ty + = 0, the equality 7V+ = N~ doesn't necessary hold for 
zero modes since D may be itself equal to zero; so Nq — Nq may be different from 
zero. Moreover, using Q7,4j), it is not difficult to see that 



Tr ( T 3 e -mA = Tr ( e -W+D^_ Tr ( e -iD-D + ) _ (7J) 



Expanding the trace in terms of zero modes and non zero modes the above relation can be 
also put in the form 

/ 



Ind (V) 



£ e"< 0) - £ e~^- 0) | + 
y s (_=o) £ (_=o) 



/ 

J2 e - tE + 0) - e~ tE( - 0) 
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leading afterwards to 

Tr (tV^ 2 ) = iV + - N . (7.8) 



(2) To get the right hand side of eq( |7.2j ), we use two more features of the index; first 
the remarkable independence of Tr (T 3 e~ tv2 ) on the spectral parameter t as explicitly 



exhibited by (|7,8|) . Second use the so called heat expansion method [45| 

Ind{V) = ^- t Y,^b n {V) (7.9) 

where the b n (£>)'s are expansion coefficients. Since this expression should be t-invariant; 
it follows that Ind(T>) = ; which by using eq(??) leads to 



b 2 (V) = Tr 



J = 2 j BdS (7.10) 



with [cr^, a v \ = 2ie^ u3 a 3 and ^- J BdS = Q producing the total magnetic monopole charge 
in discrete values inside the surface. 

7.1.2 2D fermions on lattice 

To get the topological index for 2D graphene, one starts from the analysis given in previous 
subsection; then works out the extension fermions on honeycomb. The index on a square 
lattice has been studied in |46|] ; and the one on honeycomb has been considered recently 
in [47| . Below we give a brief description of the main lines of the two constructions. 

fermions on a square lattice 
In the case of a finite L x L' square lattice where L = Na, L' = N'a with a the spacing 
lattice parameter and individual S square = a 2 , the gauge configuration is chosen as follows: 

Ax(x,y) = -uy, A 2 {x,y) = 0, F 12 =B (7.11) 

with flux 

Q =^fdxdyF 12 = ^ . (7.12) 
which, for later use, we write it as follows, 

Q = ^S B = £jQ, S = LL>, (7.13) 

We also have the following boundary conditions 

A!\y = o = a^l, + m-^n- 1 
n( x ) =e iBL ' x , o (x) u =0 = n (x) \y= L , 



(7.14) 



leading to BLL' = 2im with n integer. Notice that the discontinuity in the vector potential 
is captured by a gauge transformation. Notice also that the periodicity condition on the 
transition function VL (x) requires J^S = Q G Z; it teaches us that the field strength B and 
the topological charge Q are quantized. By following [46], the general solution of the Dirac 
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equation of the 2-dimensional fermions on the square lattice satisfying the above boundary 
conditions is given by 



kez 



(7.15) 



where H n (£) are Hermite polynomials of order n; and where the coefficients are con- 
strained by the recurrent relations Cj^ = C^_q] which result from boundary conditions and 
showing that only C , ...,C,gi 1 which can be chosen arbitrary. The eigenvalues associated 
with the ^ (x,y) } s are given by 

{E+f = 2(n + l)\B\-B , 

(E-f = 2{n + l)\B\ + B . 1 ' ' 

The link between chirality and zero modes depends on the sign of the topological charge 
Q or the sense of the external magnetic field B. For the case B > 0, the above relations 
read as (E+) = 2nB and (E+) = 2 (n + 1) B. So only the wave function with positive 
chirality has zero modes with degree of degeneracy equal to \Q\. For the case B < 0, the 
zero modes have a negative chirality with \Q\ degeneracy. 

For the case of J\f continuum flavors described by the Dirac operator in the magnetic 
background field as specified above, the index theorem reads therefore as follows 

Ind(D square )=N\Q\. (7.17) 

As such, the index of the Dirac operators for minimally doubled fermions; in particular for 
BC fermions is 2 \Q\. For naive fermions it is equal to 16 \Q\. 
honeycomb fermions 

On the honeycomb, the situation is quite similar; the main difference comes form the 
crystallographic structure. The honeycomb is given by the superposition of two sublattices 
Agra and B gra . An interesting way to parameterize these sublattices is in terms of the two 
simple roots cx\ and o: 2 of SU (3) and the weight vectors Ai, A2, A3 of its fundamental 



representations [47, 11, |49|| . More precisely sites r n in A and r' n in B are expanded as 
follows 



A: = nxoti + n 2 o.2 , n = (ni,n 2 )eZ i 
B : r' = r„ + s , d = a^- 



(7.18) 

2 

where ai, ct2 and their sum 0:3 = a± + a.2 which is also a root but not simple; can be 
taken like 

= , « 2 = (0, v / 2) , a 3 = (^,f) , (7.19) 

satisfying the usual root's properties of SU (3); in particular otf = 2. Moreover, the shift 
vector s is related to one of the weight vectors like A = s^; it can be any one of the three 
following 

si = a(l,0) , s 2 = a(-i,^) , S3 = a(-i,-^) . (7.20) 

Notice that the sum si + s 2 + S3 = 0, its captures the traceless property of the fundamental 
representation of SU(3); we also have the relations 

(si - s 2 ) = J^oli, (s 2 - s 3 ) = J| Q 2, (S3 -si) = v/|oto • (7-21) 



- 28 - 



with do = — Q3. Plaquettes in the 2D honeycomb are hexagonal with area Shexa = fl2 ^i 
so the magnetic field is quantized as follows 

B = ^—=Q. (7.22) 

In [^1 , the two lattice axes of the honeycomb were chosen as generated by the roots 
(oto, -Qi), i.e: r = xqcxq — xiQo! and the boundary conditions to have a finite translation- 
ally invariant 2D graphene lattice Lq x L\ were taken like 

F(r + L «o) =F(r) , 

F(r + L iai ) =F(r) . [ " ' 

Moreover, the link field configuration depending on these boundary conditions are given 
by 

U(r, Sl ) = e- iM ^*° , (?24) 

E/(r,s 2 ) = l , C/(r,s 2 ) = l , 

for all cells except those of the last row with xq = Lq — 1 where it is required moreover 



U(x = L - l,Xi,Sn 



i^^Loxx _ (7.25) 



Like in the square lattice, the topological index is also given by Ind(D gra ) = 2\Q\; the 
main difference is that it is given by *+E 3 ^ where E 1 , S 2 , S 3 are the generators of the 



SU (2) flavor symmetry rotating the two Dirac points; for an explicit analysis see [47]. 
7.2 Index in 4D lattice QCD 

In 4-dimensions, the determination of the index of the Dirac operator of fermions living on 
4D spaces follows the same approach as in 2D case. Let us describe briefly the main lines 
of the method. 

7.2.1 Fermions on 4D space 

To get the topological index in a SU (N) background gauge field configuration one has to 
compute both sides of the 4D analogue of eq( [7.2[ ) namely 

jV+ - N Q ~ = 2 |Q| , (7.26) 

where Nq are the number of chiral zero modes and Q the flux. Let us compute the 

topological charge Q that gives the right hand side of this relation. It is generally given 
by: 

Q = ^TT / d 4 xe^ pa Tr {J*"J»°) , (7.27) 



32vr2 JSa 

with field strength F pu = d^A u — d u A p — i [A^, A v \ valued in the adjoint representation 
of SU (N); i.e = Yla^a^u- Making a simple choice of this gauge invariant field as 
follows 

/ B x \ 
-Si 
^ ~ B 2 
\ -B 2 J 



(7.28) 
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where T is one of the generators of SU (N) normalized to Tr(T 2 ) = 2, one can determine, 
up to a gauge transformation, the corresponding gauge potentials. Notice that a gauge 
configuration that is appropriate to lattice computations is given by, 



A 2 (x,y,z,r) = 
A 4 (x,y,z,r) = 



(7.29) 



A x (x, y, z, r) = -BiyT 
As(x,y,z,T) = -B 2 zT 

Putting ( |7.28|) back into (|7.27f ) by using, 

/ -B 2 \ 
_ B 2 
-Bi 
\ B x / 

which is (anti) self dual for B 2 = =t-E>i, one gets the explicit expression of the topological 
charge Q in terms of the fields B\ and B 2 and the volume of the 4-dimensional compact 
hyper surface. We have, 

16B!B 2 



(7.30) 



Q 



32vr 2 



-Vol (S 4 ) G 



(7.31) 



7.2.2 4D hypercube 

This relation can be given a more explicit form by working by considering fermions on 
lattices. In the case of the J^D hypercube L\ x L 2 x L3 x L4, the previous relations reads 
as 



Q 



16B1B2 

32vr^ 



Vi=l 



Moreover use the fact that the field strengths I?i and B 2 are quantized as 



_ 2nN! 

we end with the following topological charge 

Q = 2N 1 N 2 , 



Bo 



2-kNo 



N h2 G Z, 



(7.32) 



(7.33) 



(7.34) 



which is independent from the nature of the 4D lattice; and is then also valid for the 4D 
hyperdiamond of J^D lattice QCD. Regarding the left hand side of ( 7.26 ); it is determined 
by solving the Dirac equation. We find 



(x,y) 



Cf , expz 



2irk x 2irk 
x H ^—z 



Li 



L, 



(7.35) 



with 



H n ,m — H n (£) x H m (£) 

e =vW\(y + t L * < 



(7.36) 



where H n (£) are Hermite polynomials; and where the coefficients Ct~ , are constrained by 



the recurrent relations 



C 



± 



&x fez Qx fez Qz ' 

showing that the degeneracy of the chiral zero mode is \Q x Qz 



(7.37) 
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8. Appendix B: Index, spectral flow and filling factor 



8.1 Index theorem and spectral flow 

In this section, we describe the main lines of the spectral flow method to get the topological 
index which can be symbolically stated as 

Ind(D) = -spf {H) (8.1) 

where spf (H) stands for spectral flow of T~L which, reads in terms of the 4-dimensional 
Dirac operator T> = 7^ (d^ — iAu), as follows 

U (m) = 7 5 (V - m) , m E R . (8.2) 

This approach was first considered in for staggered lattice fermions where the would- 
be chiral zero-modes has been identified away from the continuum limit. Then, it has 
extended in |8], || to minimally doubled fermions as well as the naive ones by using the 
point splitting method for implementing flavored mass terms. The spectral flow method 
detects exactly the index of the would-be zero modes fixing the gauge field topology; it 
has been explicitly checked numerically for 2D and 4D staggered, minimally doubled and 
naive fermions. Below, we will mainly focus on 4D; but the results can be extended to all 
even dimensions. 

The key idea of the spectral flow method relies on defining a hermitian spectral hamiltonian 
depending on two basic things. (1) The Dirac equation for zero modes namely 



DM / <f> 
-D J I X 



(8.3) 



where each block is a 2 x 2 matrix and where the operators D and are as in 
and 

*-(*\ (8.4) 
being the Dirac spinor in 4-dimensions with the two chiral components ^f± given below, 

V + =(A, tf_=(°Y 75 * ± = ±* ± . (8.5) 

(2) the spectral hamiltonian operator T~L (m) = 7 5 (T> — m) depending on a real spectral 
parameter m which can be thought of as mass. In matrix notation, we have, 

- -(v:) . « - (_**) ■ 

Notice also the traceless property of the spectral hamiltonian operator 

TrU = 0, (8.7) 
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which, a priori, should be independent of basis change in the spinorial representation space, 
turns out to play a crucial role since the sum of its eigenvalues should be equal to zero. 
Notice also that for the Dirac zero modes \E', we have 7i^> = -m7 5 f and so a zero-mode 
of D with ± chirality is also an eigenmode of % with eigenvalue A (m) = ±m as given here 
below 

n^± = =fm^>±. (8.8) 

These eigenvalues cross the axis A (m) = with two possible slopes ±1 at m = 0. Moreover, 
from the property % 2 = D^D + m, it follows that that A (m) = ±m are the only eigenvalues 
of % that cross the origin at any value of m. Therefore the spectral flow of % (m), defined 
as the net number of A(m)'s that cross the origin, counted with sign ± depending on 
the slope of the crossing, comes entirely from eigenvalue crossings at m = and equals 
N- — N + = —Ind (T>). Numerical results showed that the topological index is indeed given 
by minus the spectral flow hamiltonian; for more details see || |34| . 

8.2 Filling factor and chiral anomaly 

First recall the expression of the filling factor v in terms of the number Nf of fermions and 
the flux number N^, 

N f / s 

In 2-dimensions, we found that the filling factor u 2 D = 92D x (2-/V + 1) where g2D is some 
degeneracy factor giving the number of species and their quantum numbers. The integer 
(2iV + 1) is the sum of two contributions (N + |) coming from electrons and \N + ^) from 
holes. The two half integers | are associated with the fundamental state N = where the 
conducting (electrons) and valence (holes) bands touch; this state is a chiral zero mode as 
shown on eq (|2.19|) . Notice that the background filed is given by 

with flux $ = BS through a surface S. The same situation happens in 4-dimensions where 
we have found that v^y) = 54D x (2A r i + 1) (2iV2 + 1). The main difference is that in 4D 
the background field involves two kinds of magnetic fields 



T {AD) 



( Bi \ 
-Bi 
B 2 
\ -B 2 J 



.11) 



and then two fluxes $1 = B\S and $2 = B 2 S. The numbers (2N± + 1) and (2A^ + 1) are 
then associated with the fluxes $1 and $2- Here also the zero modes are chiral and so 
contribute to the topological index. 

More explicit relations can be written down by working on lattices on which the magnetic 
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fields and the fluxes are quantized as follows, 



Lattice 


magnetic fields 


fluxes 


2D 


i _ LiL 2 

B ~ 2-kN 1 


Q2D = 


N x 


4D 


1 B\L 2 1 L3L4 

B 1 1 B 2 2-kN 2 


Qad = 


2N 1 N 2 



3.12) 



On the other hand, using the fact that the momenta j> M = Hku of a particle of coordinate 
(x^) in background fields is given by the gauge covariant derivatives iD^ = id^ + \T ilv 



x 



we learn that in the large magnetic field this momenta (wave vector) is dominated by the 
cyclotronic term ^F^ v x v . This limit leads to 



\x»x v ] 



.13) 



where Q^ v is as in eq( [4.4| ). This relation is a typical phase space relations; it teaches us that 
in presence of a strong background field the space gets discretized into fundamental cells 
of area l 2 B = AQ^. In 2-dimensions, we have, up to a normalization factor, the following 
non commutative geometry relation 



[x,y] 



1 

B 



UL 



1^2 



2irN 



3.14) 



where the second equality, which is valid for N 7^ in agreement with the large B limit, 
follows from ( 8.12j ). From this relation, we learn the area of the fundamental cell 



I 2 
L B 



UL 



1^2 



2ttN ' 



3.15) 



and so one is left with N electrons (N holes) coupled to the quantum flux leading to the 
filling v* 2D = g2D x 2N. But this is not exactly the computed value V2D = 92 D x (2-/V + 1); 
the latter may be then viewed as the quantum version of v\ D and, due to the chiral anomaly, 
corresponds to shifting the magnetic length as 



UL 



1^2 



B 



2tt {N+\y 



3.16) 



which in turns corresponds to shifting the topological charge as Q = N N + ^. 

In 4-dimensions, the situation is quite similar to 2D case; except that here we have two 

kinds of commutations relations 



[x,y] 
krl 



B 1 

i 

B 2 



and so two kinds of magnetic lengths namely 



L\L 2 
2-KNi 1 
L3L4 
2ttN 2 ' 

LiL 2 

2ttN! 



and Z| 2 



L3L4 

2nN 2 



(8.17) 



leading after- 



wards to a classical relation of the filling factor v\ D = ANiN%. However to get the 

right expression of the filling factor namely v^d = 54D X 4 [Ni + |) (^2 + \), one has to 
shift the topological charges and the magnetic lengths as in ( |S.16D ; this behavior is also a 
manifestation of the 4D chiral anomaly. 
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